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Abstract. We have studied using the KLUN sample of 
5171 spiral galaxies having Tully-Fisher distance moduli, 
the average radial space distribution of galaxies out to a 
distance of about 200 Mpc (for H = 50 km s" 1 Mpc" 1 ). 
One motivation came from the debate on the fractal di- 
mension T> and maximum fractality scale A max of the 
large-scale galaxy distribution (Davis 1997, Guzzo 1997, 
Pietronero et al. 1997). A specific recent study is the 
3-dimensional correlation analysis of the all-sky LEDA 
data base by Di Nella et al. (1996) who concluded that 
the galaxy distribution is fractal up to scales of at least 
300 Mpc, with fractal dimension 2.2. One would expect 
to see a signal of this result in the radial space distribution 
of the all-sky KLUN sample. We have studied this question 
with a new method based on photometric TF distances, 
independent of redshift, to construct the number density 
distribution. 

Our main results are: 

1. While scattered below about 20 Mpc, at larger dis- 
tances the radial distribution starts to follow, in 
terms of distance modulus /itf, the law log A = 
(0.46±0.01)/i+const., using diameter TF relation, and 
log A = (0.40 ± 0.01)/i+const. for magnitudes. These 
are the predictions based on fractal dimensions 2.3 
and 2.0, respectively. These radial density gradients 
are valid up to the limits of KLUN, or about 200 Mpc. 

2. We have tried to understand the derived radial den- 
sity behaviour as a result of some bias in KLUN or our 



analysis, however, without success. Numerical simula- 
tions have shown that the method itself works, though 
it somewhat underestimates the radial distribution ex- 
ponent. If the density law is caused by incompleteness 
in the diameter limited KLUN sample, then the in- 
completeness should start at widely different angular 
diameters c?25 for different values of rotation parame- 
ter logVM; which would be quite unexpected. On the 
other hand, if the derived distribution is correct, the 
completeness is good down to c?25 = l'-6, as originally 
intended and previously concluded. 
3. If correlation studies favoring long scale fractality (200 
Mpc or more) and T> w 2 are correct, the position of 
our Galaxy would be close to average in the Universe, 
with the galaxy density decreasing around us according 
to the expected law (Mandelbrot 1982). 

Key words: Galaxies: spiral - Galaxies: distances and 
redshifts - Cosmology: large-scale structure of Universe 
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1. Introduction 

The present work is motivated by the longstanding prob- 
lem in observational cosmology, concerning the structure 
of galaxy distribution: is the Universe homogeneous or 
fractal on large scales, and what is the scale where ho- 
mogeneity begins? Paturel et al. (1994) and Witasse & 
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Paturel (1997) left open the possibility that the appar- 
ent partial incompleteness in some galaxy samples could 
actually reflect inhomogeneity or fractal distribution. 

The completeness of the KLUN-sample of spiral galax- 
ies was studied by Paturel et al. (1994) from the dis- 
tribution of the angular sizes (logd 2 5, with isophotal 25 
mag/arcsec 2 diameters c?25 expressed in units of O'.l). 
Separating the sample into parts with and without the 
plane of the Local Supercluster, it was concluded (on 
the assumption of homogeneity outside of the LS) that 
the sample is almost complete down to logefos = 1-2 (or 
^25 = 1'.6). In fact, it is a general problem, how to deter- 
mine the completeness of a sample, if there is no a priori 
quarantee that the galaxy counts in a complete sample 
should follow, say, the distribution predicted by a homo- 
geneous distribution of galaxies. 

We continue here the study of incompleteness (the 
knowledge of which is crucial for our applications of 
KLUN), and attempt simultaneously to determine the all- 
sky average radial distribution of galaxies. This is done 
by a method which utilizes the TF distance moduli avail- 
able for every KLUN galaxy. As the question of the radial 
distribution is closely connected with the problem of the 
scale and dimension of fractal distribution, we give below 
a short resume of the state of this field. 

The early history of determination of the fractal di- 
mension and the maximum scale of fractality is described 
by Peebles (1980): the main result was that T> w 1.2 and 
A mal < 20 Mpc 1 . The first evidence for a universal fractal 
distribution of galaxies with fractal dimension D«2 was 
given by Baryshev (1981) from number counts and virial 
mass distribution arguments. The same T> « 2 inside the 
Virgo supercluster was indicated by Klypin et al. (1989) 
from analysis of fractal properties of superclusters, with 
A mal < 40 Mpc. Later in several works a double power 
law behaviour was found so that at distances r < 6 Mpc 
V « 1.2, but for 6 Mpc < r < 60 Mpc there is V w 
2.2 (Einasto 1991, Guzzo et al. 1991, Calzetti et al. 1992, 
Einasto 1992, Guzzo 1997). 

In recent years several extensive galaxy redshift sur- 
veys have become available (CfAl, CfA2, SSRS1, SSRS2, 
Perseus-Pisces, LEDA, LCRS, IRAS, etc., see the review 
by Sylos Labini et al. 1998a). Detailed analyses of the sur- 
veys have been made by different quantitative statistical 
methods which have revealed the fractal structure of the 3- 
dimensional galaxy distribution at the scales correspond- 
ing to the surveys. The fractality implies that around 
any galaxy the number density behaves with radius r as 
n(r) = kr^-V for r < \ max . 

However, considerable debate has arisen on the value of 
fractal dimension T> and the maximum correlation length 
Amaa; (see e.g. Guzzo 1997, Sylos Labini et al. 1998b). 
One may even speak about two schools. According to the 

1 All distances in this paper are based on Ho = 
50 km s" 1 Mpc" 1 . 



first one (see Davis 1997), using the standard correlation 
function analysis, fractal structure exists within the scale 
interval 0.02 Mpc < r < 20 Mpc with dimension V w 1.2, 
while on the scales larger than 40 Mpc the Universe be- 
comes homogeneous. According to the second school (see 
Pictronero et al. 1997, Sylos Labini et al. 1998a), using a 
more general statistical approach which includes the cor- 
relation function method as a particular case, the fractal 
structure exists at least up to 200 Mpc and has dimension 
T> « 2. The debate shows the need for complementary 
approaches using other data and methods. 

Recently, Di Nclla et al. (1996) analyzed the correla- 
tion properties of the galaxy space distribution using the 
LEDA all-sky data base containing over 36000 galaxies 
with known redshifts. They concluded that the distribu- 
tion has a fractal character up to the scales of 300 Mpc 
with the fractal dimension T> « 2.2. In a recent discussion, 
Guzzo (1997) concluded that the galaxies are clustered in 
a fractal way on small and intermediate scales, and the 
turnover to homogeneity may occur somewhere between 
200 and 400 Mpc. Taken at face value, such results pre- 
dict that in an all-sky sample like KLUN, one should see 
a radial galaxy number density behaviour which deviates 
from the uniform distribution. In terms of the distance 
modulus uniformity corresponds to the differential law 
log A = 0.6/x + constant for the number N of galaxies in 
the distance modulus range fi ± l/2/i, when \i has been 
derived by the "standard candle" method. In a distribu- 
tion with fractal dimension T>, the expected density law 
around any galaxy is log A = 0.2V^i + constant, or with 
V = 2, we have the "0.4/x-law". 

There are a few noteworthy aspects in the present 
study. We concentrate on investigating how the radial 
density behaves around us, while the analyses mentioned 
above are usually based on the average behaviour around 
the different galaxies in the sample, which reveals the ac- 
tual fractal law. However, there is special importance in 
seeing directly the behaviour around our individual van- 
tage point. One thing is that we are naturally eager to 
see and check in concrete manner the "strange" predic- 
tion of fractality, or the thinning of average density, if 
our Galaxy occupies an average position in the fractal 
structure. On the other hand, Sylos Labini et al. (1996) 
have pointed out that in order to see the fractal behaviour 
around one point, one needs an all-sky sample with a suf- 
ficiently faint magnitude limit (m > 14) in order to avoid 
so-called finite- volume effect inside a few Voronoi-lengths. 
The KLUN sample can thus be used as a kind of first firm 
step for deriving the radial density law, which may be 
later straightforwardly extended deeper when the sample 
is gradually increased. 

In order to derive the density law, we use a new method 
based on photometric (TF) distances. It has two special 
advantages. First, we can use almost the whole sample for 
constructing the density law. Second, unlike the methods 
based on redshift distances, we do not require a local veloc- 
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ity field model. In view of the recent suggestion (Praton et 
al. 1997) that working in the rcdshift space may enhance 
or even produce structures of the "Great Wall" type (and 
hence influence the correlation studies and presumably, 
the derived value of V), it is especially important to use 
photometric distances in complementary studies. 

2. The sample and the method 

2.1. KLUN 

The KLUN galaxy sample and the quality of the measured 
quantities of its members have been recently described by 
Theureau et al. (1997b). Originally intended to be com- 
plete down to the isophotal diameter 1'.6, it now consists 
of 5171 galaxies covering the type range Sa - Sdm, hav- 
ing redshifts and the measured 21cm-line widths for the 
Tully-Fisher relation. 4577 galaxies have also the Bt mag- 
nitude. 

For calculation of distance moduli fi, we use the direct 
TF relations established by Theureau et al. (1997b) and 
the same corrections for type effect, internal extinction, 
and galactic absorption. For example, the diameter TF 
relation gives the distance modulus \i as 

[i = 5(log£>25 - logd 2 ° 5 + const.), 

where logics = a log Vyi + b{T). In this formula, d® 5 is the 
angular diameter corrected for galactic absorption and in- 
ternal extinction, Vm is the maximum rotational velocity 
measured from the 21cm-line width of the galaxy, b(T) is 
the TF zero point for galaxy type T, and a is the slope of 
the TF relation. 

2.2. The method of TF distance moduli 

We apply here a method outlined in Sect. 3 of Teeriko- 
rpi (1993), where the TF distance moduli can be used, in 
principle, up to the limits of the sample, in order to derive 
the form of the density law. The method utilizes galaxies 
in narrow logVM ranges to construct several differential 
distributions of fi, one for each log Vm range. These distri- 
butions are displaced relative to each other in ^t, because 
of different log Vm (absolute magnitude) and the common 
diameter limit (or more exactly, the common incomplete- 
ness law). They have also generally different amplitudes, 
because different logVM ranges are populated by differ- 
ent numbers of galaxies, though their forms are similar if 
inside the ranges the galaxies are evenly distributed. 

One expects that above some logc?25, the sample is 
complete, hence below certain limits Atii m ((log Vm)) each 
distribution reflects the true average space distribution. 
We illustrate this with Fig. 1 (see also Fig. 1 of Teeriko- 
rpi 1993), which helps one to understand why division 
of the sample into logVM-ranges allows one to separate 
the effect of true space density distribution from the ef- 
fect of incompleteness in the distribution of fi. After hav- 
ing identified this grange for the subsample with smallest 



a. before normalization 
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b. after normalization 




Fig. 1. Schematic illustration of how the distribution of TF 
distance moduli /itf differ for two narrow log Vm ranges sepa- 
rated from each other by A log Vm- Here it is assumed that the 
true space distribution produces a constant slope in the log N 
vs. /jtf diagram (corresponding to density law p(r) oc r~ a ). 
Small-log Vm subsample reveals the correct density law up 
to fi^ F where the incompleteness starting at logd 1 ^™ begins 
to affect. Large-log Vm subsample shows the same behaviour, 
though shifted towards larger /itf by Api = 5a A log Vm- In 
the present method such distributions are in addition shifted 
vertically so that the complete parts come together forming 
a common envelope which follows the true space distribution 
law. Points A and B correspond to the same angular size d' 2 s 
in the incomplete region. After normalization they also define 
the correct density gradient. 



(log Vm), hence smallest /xu m , one can normalize the other 
subsamples to this one, shifting their log N vs. \x distribu- 
tions vertically, so that up to ^n m all have the same av- 
erage normalized number counts. This is allowed because 
for the subsamples with larger (logVivi), the complete \i- 
ranges extend beyond the /in m corresponding to smaller 
(log Vm)- 




Fig. 2. logVM n distribution of the KLUN sample. Bins are Fig. 3. Differential logiV vs. /itf (diameter) distributions for 
separated with dotted lines. The dashed lines show the average different log Vm° ranges (see the text), 
values in each bin. 



The common envelope of thus normalized distributions 
should follow the space distribution law up to the ^r lm for 
the subsample with the largest (log^M)- An interesting 
feature of the method is that one may construct the space 
distribution still farther, as is shown by Fig. lb where the 
points A and B corresponding to the same value of log d' 25 
(< log<i 2 ™) in the incomplete region are connected by the 
line having the correct space density slope. The method 
is also quite simple in essence and allows one quickly to 
"see" the run of the density law and the factors influencing 
its derivation. Of course, working with real data is more 
problematic than the ideal theoretical case. For instance, 
one cannot use indefinitely small log Vm ranges. 

2.3. Special aspects 

Actually, one must use the normalized log Vm when mak- 
ing the described division of the sample: 

log V M n - log Vm - (logrf 2 °, - logd 25 )/a + (b(T) - b(6))/a 

The concept of normalization has been introduced and dis- 
cussed in a whole series of studies connected with the so- 
called method of normalized distances (see e.g. Theureau 
et al. 1997b). In this manner, galaxies of different types 
and with different inclination and galactic absorption cor- 
rections, but with equal normalized logVM 11 , have the 
same effective limiting diameters. The /XTF-distributions 
become thus better defined, and more closely follow the 
ideal distribution determined by the space density and in- 
completeness curves. In other words, for a given logVM 11 , 
the /xtf scale reflects faithfully the log cfes scale. 



Finally, we note that the method is in principle insen- 
sitive to the "patchy incompleteness" (Guzzo 1997) that 
may be present in such somewhat non-uniform collections 
of data as LEDA (the host database of KLUN which, how- 
ever, was intended to be a homogeneous part of LEDA): 
In different regions of the sky the limiting diameters (mag- 
nitudes) may be different. Let us consider two regions, 1 
and 2, where incompleteness starts at logd^ and logc£ 2 5j 
respectively, and suppose we cannot identify these regions 
from the data. The summed relation, log (Ai + A 2 ) vs. 
for a constant log Vm still follows the correct slope at small 
distances (Aq + N 2 = Aq + c ■ Aq = (1 + c)Aq, where the 
constant c depends only on the ratios of limiting angular 
sizes and region areas) while now the incomplete part of 
the curve is deformed. However, for different logVM the 
form of the curve is similar and the situation is analogous 
to what is depicted by Fig. 1. 

3. Search for the density law 

3.1. Diameter TF distance moduli 

We show the procedure in a transparent manner which 
allows one to see the steps taken and to recognize the 
impact of possible systematic errors. 

First, we have divided the sample into five logVM 
ranges according to the normalized log Vm": 1-5 - 1.7 - 1.9 
- 2.1 - 2.3 - 2.5. These intervals cover practically all the 
sample and the middle one is centered on the maximum of 
the log^M distribution (Fig. 2). Inside the extreme ranges 
the distributions are not symmetric around the mean, due 
to the systematic effect described in Sect. 3.3. 



Teerikorpi et al.: The radial space distribution of KLUN-galaxics 



5 



3. 

O 




Fig. 4. As in Fig. 3, but now the distributions are arbitrarily shifted on the /iTF-axis. 



The distributions of /itf in the log Vm" intervals are 
shown in Fig. 3. The bins on /ztf should have the width 
5 • 1.082 • 0.2/fc, (k =1,2,3,...) in order to collect the same 
l°g^25 ranges on the individual /xtf distributions as ex- 
pected from the logVM intervals (0.2) and the slope of 
the diameter TF relation (1.082). Selecting k = 2, we get 
suitable bin size A/z = 0.541. For clarity, Fig. 4 shows the 
distributions artificially shifted on the /LtTF-axis, revealing 
their similar forms. 

The next step is the normalization of the numbers 
logA^. For this we use the following procedure: 

1. Select the complete part of the logc?25 distribution. 
The putative limit is logons = 1-2, as expected from 
previous studies. 

2. Find the corresponding part of the first subsample 
(with smallest logVM n ) and count the total number 
of galaxies below and at /iii m (l) = 30.5: N(l, /Kii m (l)). 

3. For the second subsample calculate N(2, /zii m (l)). 

4. Shift curve 2 by AlogiV = log iV(l, Miim(l)) - 
log7V(2, Mlim (l)). 

5. Normalize the remaining curves using the same proce- 
dure with increasing /ii; m by the shift between the sub- 
sequent curves. In this manner we utilize progressively 
deeper complete parts. Then the cumulative correction 
for the j:th curve is 

i 

A j \ogN = Y J [logiV(i-l,Atiim(i-l)) 

8 = 2 

-logAT(i,/ii im (i - 1))] , 

where /zn m (i) = /in m (i — 1) + 0.541. By restricting 
the steps in this manner, we keep progressively farther 
from the putative limit logd,25 = 1-2, while having an 



increasing number of galaxies available for normaliza- 
tion. 

Figure 5 shows the resulting normalized composite dis- 
tribution. One may note two particular things: There is 
first a rather scattered run of points up to about /xtf = 31 
which can be roughly approximated by a 0.6/xtf line. Af- 
ter that the common envelope assumes a shallower slope 
of about 0.46. 

Using now the points above log (I25 = 1.2 for each curve 
and requiring that each point contains more than 12 galax- 
ies, by which the noisy //tf ranges are removed, we draw a 
straight line through them. For this line, there are points 
up to /xtf = 34. We use this line, derived by simple least- 
squares technique, as a reference for further study of the 
density law using information from the incomplete parts of 
the /i - distributions. This "envelope line" is shown in Fig. 
6 together with points from the incomplete parts where we 
again use only those which contain more than 12 galax- 
ies The average (logc^) is written besides each symbol. 
One can easily discern the sequences corresponding to the 
same positions in the logd 2 5 distributions (c.f. points A 
and B in Fig. 1). These averages do not vary much along 
the same sequence. Comparing with the line of slope 0.46 
inserted, one sees that the sequences follow rather well 
this slope, similarly as was found above for the points in 
the complete part (especially above /xtf ~ 31). 

We tested the completeness limit by assuming limits 
at logc?25 = 1-3 and 1.4 instead of the above logc?25 = 1-2- 
With these more conservative values we can be more cer- 
tain of the completeness of the sample, however at the 
price of reduced number of data points. We got slopes 
0.46 (±0.01) and 0.43 (±0.02) for limits 1.3 and 1.4 re- 



6 



Teerikorpi et al.: The radial space distribution of KLUN-galaxics 




/2rj.p 

Fig. 5. Normalized composite distribution, constructed from 
the distributions of Fig. 3 by the steps described in the text. 
Note the appearance of a good envelope line. The inserted dot- 
ted lines have slopes 0.6 and 0.46. logiV = levels are conve- 
niently identified as horizontal parts. 
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Fig. 6. The envelope line from Fig. 5 (shown as dotted line), 
together with points from the incomplete parts. The dots repre- 
sent the envelope points and other symbols refer to incomplete 
parts. The numbers besides the latter give (logffes) and allow 
one to easily recognize the sequences following a slope close to 
0.46 (envelope line). 




Fig. 7. The line of slope 0.6 forced to go through the first 
normalization point. In this case the incompleteness is seen to 
increase together with log Vm ■ 

spectively, but the "envelope lines" in both cases were not 
as well defined as before. For now, we are more satisfied 
with the completeness limit at logefos = 1.2 giving us the 
slope 0.46 (see also Sect. 3.4 for another test, using simu- 
lated galaxy samples). 

Naturally, all this is interestingly close to 0.44 pre- 
dicted on the basis of fractal dimension T> w 2.2, as ob- 
tained e.g in the Di Nella et al. (1996) correlation analysis 
of LED A. However, we are well aware that systematic ef- 
fects may be involved, when we are working rather close 
to the completeness limit of the KLUN-sample. We have 
attempted to find any systematic error that could explain 
the shallower than 0.6 slope and have also looked what 
happens when one forces a line of slope 0.6 to start at 
A*tf = 31. In Fig. 7 this has been done, and the line which 
rather well describes the (scattered) run of data below 
Mtf =31, systematically deviates from the envelopes of 
the normalized distributions at larger /iTF • This behaviour 
implies that if the deviation from the 0.6/z-law is due to 
incompleteness problem, the incompleteness in apparent 
size starts for different logVvr ranges at different logc?25- 
This is something that we cannot understand, because the 
needed effect is quite large. For the largest log Vm it means 
that incompleteness starts around \0gd25 = 1-7, while for 
the smallest logVM such a limit would be around 1.2 (in 
terms of magnitudes this corresponds to a difference of 
about 2.5 mag). Still another way to state the problem is 
that if we try to shift the maxima of the distributions to 
follow the 0.6^i-slopc, the curves are everywhere separated, 
there is no normalization and no common envelope. A 
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clear argument against such a large effect comes also from 
the method of normalized distances used in Theureau et 
al. (1997b): in the "unbiased plateau" produced as a part 
of the method, one should readily recognize such differ- 
ences in the selection functions of galaxies with small and 
large logVM- Also, from the manner of how KLUN was 
created as a diameter limited sample, independent of any 
considerations of log Vm, there is no reason to expect so 
significant dependence of the selection on log Vm • For in- 
stance, a recent analysis of the H I line profile detection 
rates at Nangay radio telescope by Theureau et al. (1997c) 
does not give any indication that large log Vm galaxies are 
significantly underrepresented in KLUN. 

Finally, if erroneous, the present slope 0.46 is just by an 
accident very close to ones obtained by several quite dif- 
ferent correlation analyses (e.g. Sylos Labini et al. 1998a). 

3.2. Magnitude TF distance moduli 

The B-magnitudc TF relation has smaller scatter than the 
diameter relation, which makes it tempting to use also it 
in this study even though the numbers are then smaller 
and the selection properties of the resulting sample are 
complicated by the fact that KLUN has been originally 
selected on the basis of apparent size. However, if the dis- 
tribution of the galaxies of different log Vm ranges is the 
same in the log g?25 ~ Bt plane, one can use the magnitude 
relation in a similar manner. 

We do not go through the steps in such detail as for 
the diameter distance moduli. Figure 8 shows directly the 
normalized composite diagram. Because of the underlying 
diameter limit, the incompleteness begins farther from the 
maxima than in the case of diameters, and the complete 
part of the envelope is now shorter in /ztf- The envelope 
line constructed from points below the putative magni- 
tude limit B = 13.25, is also shown. It has now the slope 
0.40. Because the numbers of galaxies at each point of this 
diagram are smaller than in Fig. 5, the error of the slope 
is slightly larger than when using diameters (c = 0.012, 
compared to a = 0.010 with diameters). Inspection of the 
points in the incomplete part shows that the sequences, 
analogous to what was discussed for the diameter moduli, 
have similar slopes. However, now the achievable ^-ranges 
and numbers are smaller, and we do not show them sepa- 
rately. 

3. 3. Systematic error caused by finite log Vm ranges 

One can see a systematic effect in this method where 
we are forced to use finite log Vm ranges instead of ideal 
infinitesimals. It comes from the fact that galaxies are 
not quite similarly distributed inside the different logVM 
ranges. The distribution of log Vm has a maximum. Be- 
cause of the form of the log Vm distribution, inside the 
small - log Vm interval one expects an increasing number 
density of galaxies towards the edge with larger logVM, 




30 35 



Fig. 8. Normalized composite distribution for B-magnitude 
TF distance moduli in different log Vm" ranges. The slopes 0.6 
and 0.4 are shown as well as the envelope line. 
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Fig. 9. Normalized composite distribution for diameter TF 
distance moduli, constructed as Fig. 5, but using intervals of 
0.1 in logVM n in the range 1.7 - 2.3. Note the appearance of 
a good envelope line. The inserted dotted lines have slopes 0.6 
and 0.41. 
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while for the large - logVM range this trend is reversed 
(see Fig. 2). In Fig. 6 one sees that the averages (logefes) 
within different sequences do not vary very much which 
suggests that the effect is actually not very important. 

In order to check whether decreasing the interval size 
influences the result obtained above, we made an experi- 
ment whereby log\4i n interval was reduced to 0.1, in the 
range 1.7 - 2.3, where the numbers of galaxies remain 
large enough. Now the slope of the envelope line is 0.41, 
the furthest point of the line being slightly above 100 Mpc 
(Fig. 9). Again, the diminished number of galaxies make 
error of the slope larger, a = 0.015, while a = 0.010 in 
Fig. 5. 

3.4- Numerical experiment using simulated galaxy sample 

In order to check further the reliability of the used 
method, we have made numerical experiments with simu- 
lated galaxy samples. In making these tests, we have kept 
in mind the following points: 

1. Because individual distance moduli have considerable 
scatter and because we are interested in the all-sky av- 
eraged behaviour of radial density, the present method 
and the available data allow one to derive a quite 
smoothed-out view of the space density around the 
Galaxy. 

2. The all-sky average, as derived by the present method, 
does not make a difference between a random distribu- 
tion of galaxies with a radial density variation and a 
fractal distribution with the corresponding fractal di- 
mension. Hence, for the purpose of testing systematic 
errors in the method, it is sufficient to consider simu- 
lated distributions, where randomly scattered galaxies 
have a smooth radial density variation. 

3. The all-sky averaged radial distribution inside a frac- 
tal structure is statistically the same around all 
galaxies, which gives special motive for applying the 
present method. In non-fractal structures, such as 
supercluster-void network (e.g. J. Einasto et al. 1997), 
the radial distribution depends on the position of the 
observer (though in many cases the presence of the 
plane gives an apparent D « 2, also reflected in corre- 
lation function analysis where actually an average of all 
the observers is taken, see the models by Einasto 1992). 
It is intended to extend the present method to study 
large individual structures in specified regions of the 
sky detected previously with redshift-distances (such 
as the Great Wall). Such applications will need spe- 
cially tailored simulations which show how the method 
draws the density distribution curve, say, through a 
narrow plane of galaxies. 

In the experiments we started with large number of 
galaxies (say 10 5 ), for which we alloted radial distances 
using an input value of the radial density gradient a = a- m 
in the density law log ./V = a/j,+ const. For each galaxy we 
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Fig. 10. Test of the apparent diameter completeness limit 
logdn m using simulated galaxy samples. Dotted line shows the 
input value of density gradient a = 0.47, which corresponds 
to fractal dimension D = 2.35. For more complete samples 
(smaller log dn m ) the observed D approaches the input value. 
Error bars are mean deviations (2<r) for 1000 simulations. 
logdn m = 1.2 is the completeness limit assumed for KLUN. 



chose a random absolute diameter from a gaussian distri- 
bution with log -D a bs, o = 1-3 and a = 0.18 (-D a bs in kilo- 
parsecs) . Even though the log -D a bs distribution in reality 
hardly is gaussian, it is not too far from the truth for the 
KLUN sample. All the galaxies having apparent diameter 
larger than limit log dn m were included in the "observed 
subsample". To make this observed sample more realis- 
tic we also allowed in some galaxies below the diameter 
limit, percentage of included galaxies being progressively 
smaller further below the du m . Each of these observed 
galaxies were then given a rotational parameter logVM 
by inverse TF relation (logVM = a'logZ? a b s + b') with 
a' = 0.5, b' = 1.13 and a gaussian dispersion a = 0.15. 
These values resemble results of the recent KLUN study 
(Theureau et al. 1997a). Now the log N - /iTF graph (dis- 
tance moduli from direct TF relation) for the observed 
subsample could be investigated as above was done for 
the real KLUN sample. 

First we selected logG?n m = 1.2 and varied the density 
gradient a; n . For each a m we chose the initial number of 
galaxies so that the number of observed galaxies was the 
same as in KLUN sample. Then we calculated the slope of 
the "envelope" line in the logiV - ^tf graph getting the 
observed density gradient a bs For each a m we repeated 
the simulation 1000 times to get (a bs) with error bars. 
For every a m the resulting (a b s ) was smaller than a- m . 
However, this tendency was not large enough to explain 
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the deviation of our observed slope (in this section we use 
as a reference point a = 0.44, which is a weighed average of 
the slopes obtained with diameters and magnitudes) from 
homogeneous galaxy distribution (a = 0.6). The simula- 
tions showed that for a obs = 0.44, a in = 0.47 ±0.04 (2cr 
errors). In terms of fractal dimension (a = 0.2P) we can 
say that the observed value V = 2.2 can be affected by 
our methods so that the true value is T> w 2.35 ±0.20. 

We then tested the effect of the completeness limit 
by varying log c?ii m , while keeping the other parameters 
fixed. Figure 10 shows how T> ^ s — > T> ln when the sample 
gets more complete. Also, with smaller completeness limits 
the number of observed points increase, and the errors get 
smaller. This emphasizes the importance of expanding the 
database to make it deeper and more complete. 

4. Discussion and concluding remarks 

As it appears from the preceding analysis, we have ar- 
rived at the conclusion that KLUN galaxies seem to in- 
dicate a "thinning" of the galaxy universe from ^tf ~ 
31 up to « 36. In order to see the actual range, the 
distance moduli must be corrected for the Malmquist 
bias of the first kind which in the case of V = 2 be- 
comes V/'ix 1.382 a 2 M = 0.92cr^ for the magnitudes and 
V/3x 5 x 5 x 1.382of og£) = 23.05af og D for the diameters 
(see e.g. Teerikorpi 1997). Using a — 0.5 mag and = 0.15 
for magnitudes and (log) diameters, respectively, these cor- 
rections make the maximal distance probed by both the 
diameter and magnitude relations about 200 Mpc, corre- 
sponding to the first point beyond the maximum of the 
highest logVM distribution. 

This result, using photometric distances independently 
of redshifts, and an all-sky sample, offers a complemen- 
tary piece of evidence for the debate on the dimension 
and maximum scale of fractality. The value V w 2 which 
has appeared in several correlation studies (Klypin et al. 
1989, Einasto 1991, 1992, Guzzo et al. 1991, Calzetti et al. 
1992, Di Nella et al. 1996, Guzzo 1997, Sylos Labini et al. 
1998a), is here seen as a density law up to 200 Mpc, the 
limit of the KLUN sample. Taken alone, the present study 
does not prove any large scale fractality, it only suggests 
the existence of such an average density law around us 
(which in a restricted distance range could be caused by 
filamentary structure and voids or a large flattened sys- 
tem; Einasto 1992, Einasto et al. 1994, Paturel et al. 1994). 
On the other hand, within the framework of a large scale 
fractal distribution, our position on a peak surrounded by 
a decreasing density is as expected (Mandelbrot 1982). In 
this sense, when considered together with large scale cor- 
relation studies leading to V w 2, the present result is 
consistent with the view that the Milky Way occupies a 
typical position in the fractal structure. This agrees with 
recent studies by Karachentsev (1996) and Governato et 
al. (1997) who conclude that the Local Group is a typi- 



cal representation of other small groups, both in structure 
and environmental dynamics. 

The result of our work is also consistent with recent 
results on the supercluster-void network (M. Einasto et 
al. 1994, 1997) and arguments based on the correlation 
analysis of rich clusters of galaxies (Einasto & Gramann 
1993, J. Einasto et al. 1997). These studies claim that the 
observed transition scale to homogeneity is about 200-300 
Mpc. 

In view of local structures, fractal or not, we would not 
expect any perfect 0.6/i-law. Nevertheless, we have tried to 
understand if such a strong deviation from a homogeneous 
distribution could be due to some systematic error. How- 
ever, numerical tests have shown that the method works, 
though it tends to somewhat underestimate the density 
gradient. 

Naturally, the present study must be seen as a first 
step which uses distance moduli in the described manner 
and the method must be developed along with the in- 
creasing size of KLUN. A promising prospect is offered by 
the ongoing DENIS project (collecting 200000 galaxies by 
1999) and the FORT project improving the efficiency of 
the Nancay radio telescope (see Theureau 1997). A deeper 
sample will allow us to better recognize systematic errors 
and the real space density trend. 
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